On calculating the Berry curvature of Bloch electrons using the KKR method 
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We propose and implemented a particularly effective method for calculating the Berry curvature 
arising from adiabatic evolution of Bloch states in k space. The method exploits a unique feature of 
the Korringa-Kohn-Rostoker (KKR) approach to solve the Schrodinger or Dirac equations. Namely, 
it is based on the observation that in the KKR theory the wave vector k enters the calculation only 
via the structure constants which reflect the geometry of the lattice but not the crystal potential. For 
both the Abelian and non-Abelian Berry curvature we derive an analytic formula whose evaluation 
does not require any numerical differentiation with respect to k. We present explicit calculations 
for Al, Cu, Au, and Pt bulk crystals. 
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I. INTRODUCTION 

Over the past decade it has been realized that the 
Berry curvature f2 n (k) associated with Bloch waves in 
solids can play an important role in spin and charge 
transport by electrons.—^ Consequently, a first principle 
calculation of this qunatity was highly desirable. Impor- 
tant examples where such calculations have already been 
found useful are the Anomalous Hall Effect (AHE) 3 ~ 5 
and the Spin Hall Effect (SHE)££. Particularly insight- 
ful are those which focus on the integral of r2 ra (k) over 
the Fermi surface only. Following Haldane's suggestion, 9 
it was applied in Ref. d. 

The methodologies used in these calculations are based 
on two disctinct approaches. One is the evaluation of the 
Kubo formula for the off-diagonal elements o~ xy of the 
static conductivity^^— ^ The other one uses the first 
principles Wannier representation of the Bloch states4^ 
In what follows, we present an alternative way con- 
structed within the framework of the KKR approach. ! 1 ^ 12 

Since the most interesting problems, where the above 
curvature is relevant, concern the role of spin-orbit cou- 
pling, we want to develop our approach for a fully rel- 
ativistic description of the electronic structure. To be 
more specific we recall that for the Dirac Bloch wave of 
the conventional formi 3 - 



*«k(r) = e 4k r Wnk (r) , 



(1.1) 



the connection corresponding to the geometrical phase 
7„(k) is defined as 



A n {k) = i J M J t lk (r)V k Unk(r)c!r , 



(1.2) 



where the integral is over a unit cell of the volume u>. 
Then, the corresponding curvature is given by 



n n (k) = v k x A n (k) 



(1.3) 



In the above notation n is a band index and u„ k (r) is a 
periodic four component spinor function of r. 



Here we shall demonstrate that KKR-based band the- 
ory methods are particularly well suited for the task. Our 
central point is that the KKR matrix, whose determi- 
nant is conventionally used to find the energy bands, has 
its own well defined geometrical phases, connections and 
curvatures. As well as being easy to calculate, they are 
closely related to those defined above. The root cause 
of this convenient feature is the fact that such matri- 
ces depend parametrically on the wave vector k and the 
energy E . Therefore, the geometry of their eigenvalues 
and eigenvectors, in the k and £ space, is closely related 
to that associated with the periodic part of the Bloch 
functions^ There are three factors which make the study 
of KKR matrices computationally efficient. Firstly, by 
the standards of first principles electronic structure calcu- 
lations the ranks of KKR matrices are quite small. Typ- 
ically, one is dealing with 16x16 (Schrodinger equation) 
or 32 x 32 (Dirac equation) matrices (if we assume one 
atom per unit cell). Secondly, the crystal momentum k 
enters into the computation only through the structure 
constants and their gradients with respect to k. Further- 
more, these quantities depend only on the geometrical 
crystal structure but not the crystal potential. So, they 
are readily calculated, without taking numerical deriva- 
tives, by using the so-called screened version of the KKR 
methodi 14 ! 15 Finally, the calculations can proceed in the 
constant energy mode which is particularly efficient when 
studying Fermi surface properties. 

To demonstrate the efficiency and stability of the pro- 
posed numerical procedures, we present explicit calcula- 
tions of the Berry curvature on the Fermi surfaces of Al, 
Cu, Au, and Pt bulk crystals. In a fully relativistic the- 
ory the presence of both space and time inversion sym- 
metry forces every k state to be twofold degenerate! 16 ! 17 
As a consequence, we have to deal with the so-called non- 
Abelian Berry curvature i 18 ' 19 The corresponding formal- 
ism is derived within this paper. In order to illustrate 
the significance of such calculations for the understand- 
ing of interesting physical phenomena, we computed the 
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intrinsic contribution to the spin Hall conductivity for 
Pt and Au. We compare our results with those obtained 
by other methods.—^ Although these calculations were 
performed using the usual Fermi sea integration, our for- 
malism should be especially efficient for approaches based 
on Haldane's suggestion.— According to Haldane the cal- 
culations require the Berry curvature only on the Fermi 
surface. 

We will introduce our novel theoretical framework for 
calculating the above connection and the curvature in 
two steps. In Section II we present an alternative for 
computing the group velocity 

v n (k) = V k £ n (k) (1.4) 

of Bloch electrons, without taking the partial derivative 
with respect to k numerically. In Section III we extend 
this approach to the calculation of the Berry curvature in 
both, Abelian and non-Abelian, cases. The example com- 
putations are shown and discussed in Section IV. Note 
that in Sections I-IV we use atomic units with energy 
in Rydberg. The results for the spin Hall conductivity 
are presented in Section V. We conclude in Section VI 
and the Appendix provides a detailed derivation of the 
formulas used in the calculations. 



II. CALCULATING THE GROUP VELOCITY 
FOR BLOCH ELECTRONS 

In this section we prepare the ground for our princi- 
ple task in Section III by outlining a simple, instructive 
way of computing the group velocity v„ (k) . In addition, 
we introduce briefly the relativistic KKR formalismi 21 i 22 
The basic idea for calculating the group velocity was sug- 
gested by Shilkova and Shirokovskii in Ref. [23|. However, 
our procedure will follow a slightly different route and 
hence will be described in detail below. 

To be specific with regard to notation we use that of 
Ref. [ST]- Here we restrict our consideration to the non 
spin-polarized case that means nonmagnetic systems. To 
simplify the equations, we assume one atom per unit cell. 
Nevertheless, the generalization to a lattice with a ba- 
sis is straightforward. In addition, we use the atomic- 
sphere approximation (ASA) for the crystal potential in 
the Dirac equation^ 

Then the Bloch wave corresponding to a band n can 
be expanded around a site in the ASA sphere as 

* Tlk (r)=EC3(k)$g(£ n (k);r), f2 ±) 

where 

<M£;r) = ( 9 r { fA XQ H) (2-2) 

are the scattering solutions of the Dirac equation for the 
spherically symmetric potential at the energy £ . They 



are written in terms of the large and the small compo- 
nent, where gx(£;r) and /*{£ ;r) are the corresponding 
radial functions * 21 ' 22 Here Q = {x, ^1} and Q = {—x, /i} 
are abbreviations for the quantum numbers jt and /1 
specifying the conventional spin-angular eigenfunctions 
^g(e r ) f i^ where e r = r/r. 

When the multiple scattering ideas of Korringa^i and 
Kohn and Rostoker— are invoked, one finds that the en- 
ergy eigenvalues £ n (k) are given by those combinations 
of £ and k for which the determinant of the KKR matrix 

M QQ , {£■ k) = Gqq, (£; k)At° Q , {£) - 6 QQ , (2.3) 

is zero. Note that the screened structure constants 
Gqq,(£; k)— depend only on the crystal structure while 
the screened At-matrix describes the scattering at the lo- 
cal, self-consistent effective one-particle potential. There- 
fore, Atq(£) is a function of energy £ but not of k. This 
is the separation of crystal structure and potential men- 
tioned in the introduction. Moreover, the more sophisti- 
cated, and physically more relevant, spin-polarized ver- 
sion of the theory will retain the formal structure with 
the difference that the At-matrix will be non diagonal in 
Q. 

An efficient way of finding the zeros of the KKR deter- 
minant ||A/QQ/(£;k)|| is to solve the matrix eigenvalue 
problem 

M(£;k)C n = \ n C n (2.4) 

and to search for vanishing eigenvalues A„(£;k). It can 
be performed, either in k space at constant energy or in 
£ at fixed k. In the above notation the components of 
the matrix M and the nth eigenvector C n = {Cg(k)} are 
labeled by Q. 

By means of the expansion coefficients C n , correspond- 
ing to the band energy £„(k), we could calculate the 
group velocity evaluating 

v n (k) = J *, t lk (r)ca* nk (r)dr , ( 2 .5) 

with the relativistic velocity operator col. As was shown, 
analytically, by Shilkova and Shirokovskii 2 ^ this formula 
is equivalent to Eq. (|1.4j) . However, within the ASA 
approximation used in this paper, the expression (which 
follows from Eqs. (gUJ and (123]) ) 

v n (k) = Cl(k)c&(£)C n (k) , (2.6) 

where the elements of the vector matrix a are defined as 

(a) QQ >(£) = f& Q (£;r)&$ Q ,(£;r)dr , ( 2 .7) 

does not reproduce the results of the numerical differen- 
tiation exactly. We will comment on this problem at the 
end of the current section. 

For now we turn to the central result of Ref. [23J which 
is based on Eqs. (|1.4p . (|2.4I) and derive a similar expres- 
sion. Technically, the solution would be easier if the ma- 
trix M{£ ; k) was Hermitian. However, due to the used 
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expansion, it is not. An additional transformation, dis- 
cussed by Kohn and Rostoker— and used in our previous 
papers based on Refs. [H] and Hi], can provide a Hermi- 
tian KKR matrix. 24 However, the derivation of the Berry 
curvature described in Section III would be more com- 
plicated due to the necessary normalization of the basis 
functions. Thus, for clarity, we proceed to solve Eq. (12. 4[) 
for a non-Hcrmitian matrix M(£;k). In short, we find 
the right and left eigenvectors, C n {£; k) and D n (£; k), re- 
spectively, such that the following conditions are fulfilled: 
C^Cn = 1, D^Dn = 1. We note that D^C^ oc 8 nn i, 
CiCn* ± 0, DlD n , ^ 0, and £>tC„ ^ l. Here C n and 
D n correspond to the same eigenvalue A„. A straightfor- 
ward algebra, summarized in Appendix [5] yields 

V «W = - 1^=1 " • (2 ' 8) 

This expression, being the main result of the cur- 
rent section, is similar to the one obtained in Ref. 
for the Hermitian KKR matrix. It shows that having 
found the Bloch state energy £„(k), that provides the 
zero of the nth eigenvalue A„(£;k), one can calculate 
the velocity by evaluating the above formula. For this 
purpose, the eigenvectors C n (£;k) and D n (£; k) cor- 
responding to A„(£„(k);k) = as well as the partial 
derivative of M(£; k) with respect to k are required. 
Since in the screened KKR method dM/dk can be eval- 
uated analytically, the disadvantage of taking numerical 
derivatives of the dispersion relation £ n (k) is avoided. 
Namely, it follows from Eq. ([231) that dM(£;k)/dk = 
(dG s (£;k)/dk)At s (£). Noting that, one can use the 
short range feature of the screened real space structure 
constants G s (£; R) 1 ^ to evaluate 

dG s (£;k)/dk = iJ2Re ikR G s (£;Il) (2.91 



discussed in the literature with respect to dipole tran- 
sition matrix elementsi 25 ' 26 It was shown that the ASA 
approximation causes difficulties in evaluating the off- 
diagonal matrix elements of the relativistic velocity op- 
erator. The authors of Refs. [25| and l26l resolved the issue 
by rewriting the necessary formulas to get numerically 
more stable results. The problem in evaluating the ex- 
pectation value of ca was already discussed by Shilkova 
and Shirokovskii who solved the problem by following the 
line of arguments we have adopted here. They showed 
that this method is perfectly stable. Here we confirm 
their results for the case of non-Hermitian KKR matri- 
ces. 

Finally, we point out that the method of calculating 
the Berry curvature presented in the next section uses the 
same techniques as considered above. Therefore, similar 
improvements of accuracy and stability for the numerical 
results are expected. 
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FIG. 1: The absolute value of the Fermi velocity of Cu (in 
a.u.) obtained using three different methods: a) the numerical 
derivative of the dispersion relation Vk£n(k); b) the expec- 
tation value of the relativistic velocity operator given by Eq. 
(|2.6p : and c) implementation of Eq. (|2.8jl 



III. NEW ROUTE TO COMPUTE THE BERRY 
CURVATURE 



at each k point, separately. Consequently, the only nu- 
merical derivative to be taken, by calculating the velocity, 
is the one-dimensional derivative d\ n {£ ; k)/d£ . Fortu- 
nately, this requires only modest computational efforts. 

In concluding this section, we report in Fig. 1 a com- 
parison between Vk£ n (k) as calculated by numerical dif- 
ferentiation, by the use of Eq. (|2.6p . and by evaluating 
the formula of Eq. (12.81) . The calculations are performed 
for the electron states on the Fermi surface of Cu. Sig- 
nificantly, the results based on Eqs. (|2.6p and (|2.8[) show 
a smooth appearance over the Fermi surface, indicating 
their independence on the number of k-mesh points. By 
contrast, the numerical derivative in Eq. (|1.4p strongly 
depends on the used k mesh. The other noteworthy fea- 
tures of these results are the similarities and differences of 
the velocities obtained by Eqs. (12.61) and (12.81) . A detailed 
analysis of those shows that the numerical derivative of 
Vk£ n (k) converges to the result of Eq. (|2.8p whereas for 
the direct evaluation of the velocity operator by Eq. (12 
a maximal error of 4% remains. This effect was already 



In this section the formalism for the calculation of the 
Berry curvature within the KKR method is derived. We 
start with the conventional ( Abelian) case for A n (k) and 
0„(k) (subsection A and B, respectively). Then we ex- 
pand our consideration to a general non-Abelian case 
(subsection C). 

A. The connection for u n ^(r) via >]/ n k(r) 

Clearly, the periodic part w n k(r) of the Bloch wave 
is an eigensolution of the Schrodinger or Dirac equation 
with Hamiltonian 

H k (r) = e - lkr i?(r)e lkr (3.1) 

in which the wave vector appears as a parameter. Thus, 
the arguments leading to Eqs. (|1.2|) and (|1.3p are, by 
now, conventional. 2 However, whether the Bloch wave 
itself has a geometrical phase, connection and curvature 
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in its own right appears to be a different problem. The 
Hamiltonian for ^ n k(r) does not depend on k, and the 
wave vector enters into the discussion of Bloch waves 
only by defining the boundary conditions. Although it 
has been noted2£, it was not clarified whether a slowly 
changing boundary condition is exactly equivalent (or has 
the same holonomy) to a slowly changing parameter k in 
the theory of u„k(r). 

Another comment which concerns the above discussion 
is that k of v T / n k( r ) labels not only the energy eigenstate 
but also the eigenvalues e lk ' R of the translation operators 
Tr. Therefore, it is not entirely free to act as a parame- 
ter. By contrast, u nk ( r ) is degenerate with respect to all 
translation operators and hence its k is not obliged to la- 
bel their eigenvalues. As a consequence, they are free to 
be parameters in Hk(r). In other words, u„k is not in the 
same Hilbert space as it„k' and hence they do not need 
to be orthogonal. In contrast, u n k and w m k with m =/= n 
reside in the same Hilbert space and are orthogonal to 
each other iSI 

With these remarks in mind we note that the KKR, 
as most band-theory methods, is designed to calculate 
'J'nk(r) but not it nk (r), in addition to the energy eigen- 
value £ n (k). Nevertheless, the Bloch function in the unit 
cell u>, as given by Eq. (|2.1[) , can be used to evaluate the 
connection as 

A n (k) = i J *t k (r)V k *„k(r)dr + / ^ k (r)r* nk (r)dr . 

From the point of view of the above discussion it should 
be stressed that the integrals in the above expression are 
over a chosen unit cell only and they are not the usual 
matrix elements between Bloch states. Clearly, such ma- 
trix elements would feature integrals over all the space 
with the corresponding orthogonality. In contrast, while 
integrating over a unit cell the Bloch states are not or- 
thogonal. 

The purpose of writing A n (k) in the form of Eq. (|3.2[) is 
not to attribute it to the Bloch states, but to facilitate its 
calculation using the local expansion of Bloch states given 
by Eq. (I2.1[) . As will become apparent shortly, the two 
contributions on r.h.s. of Eq. (|3.2[) correspond to different 
aspects of the problem. Therefore, it is convenient to deal 
with them separately. For each reference we shall call the 
first term A!^(k) and the second A r n (k). 

Let us use the KKR expansion given by Eq. (|2.1|) and 
the fact that the scattering states $g(£;r) can be nor- 
malized to 1 within a unit cell. Then, a straightforward 
calculation of A5(k) yields 

^(k) = ^™(k)+^(k), (3.3) 

where (a detailed derivation is given in Appendix [B]) 

A v n (k) = iv n Ct KC n = -v n Im{CiKC n } (3.4) 

with 

(K) QQ ,(£) = 6 QQ ,J^ Q (S;r)^^ldr (3.5) 

W 



and 

A« KR (k) = iClVA = -Im{Ct V k C„} ■ (3.6) 

Here the matrix A is diagonal because the angular part 
of the KKR-basis set (Eq. (|2.2p ) does not depend on en- 
ergy. Clearly, the term given by Eq. (|3.6p is similar to 
the standard formula for the connection. It is associated 
with the eigenvalue problem of Eq. (|2.4[) in the usual 
wayi and hence can be regarded as a property of the 
KKR matrix M(£\k). This term is purely real since 
C^VkCrj is a purely imaginary quantity due to the nor- 
malization C\C n — 1. The other term, A„(k), is always 
parallel to the group velocity and purely real due to the 
antihermitian property of the matrix A. 

Turning to the second term in Eq. Q3.2p and using the 
local expansion of Eq. (|2.1|) one readily finds 

A^(k) = Ct(k)FC„(k) , (3.7) 

where the vectorial matrix r is defined as 

(F) QQ , {£) = / $ Q (£; r)r<fv (8; v)dr . (3 . 8) 

Then the full connection is given by 

A n (k) = A? KH M + ^(k) + A(k) (3.9) 

together with Eqs. d3T4j) , (|3l|, and (|3~7l) . 

The next subsection is devoted to present a method for 
calculating the curvature given by Eq. (|1.3[) within this 
framework. 



B. KKR formula for Abelian Berry curvature 

It follows from Eq. (|3.9[) that the curvature can be 
considered as a sum of the following contributions 

n n (k) = n™(k) + n'(k) + n;(k) . (3.10) 

We start with the first term of r.h.s. in the equation 
above, namely fl^ KR (k) = V k x A% KR {k). This is the 
curvature associated with the KKR eigenvalue problem 
of Eq. (|2.4j) . To deal with it, we note that 

tt™(k) = iV k Ct x V k C„ = -Im{V k Ci x V k C„} . 

This is the standard form of the Berry curvature derived 
from a matrix eigenvalue problem^. However, because the 
KKR matrix M(£; k) is not Hcrmitian, the algebra from 
here on deviates somewhat from the usual procedures.^ 
In particular, the completeness relation J2 m Cm,C} n = 1 
for M(£; k) beeing Hermitian fails in our case. Instead, 
to transform Eq. (|3.11|) into a computationally conve- 
nient form, we must use 

N - -+ N -t 

v S^EL = v S^SL = I , (3.12) 

m—1 m—1 
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where as before C m and D m are right and left eigenvec- 
tors of M(£ ; k), respectively^ Here the sum is going over 
all N eigenstates of the matrix M that has a dimension 
of N x N. Substituting Eq. (|3~T2l) into Eq. (j3TTTT) the 
KKR curvature takes the following form 



ft 



KKR 



(k) = -7m{£ 



} • (3.13) 



The next move is to eliminate the derivatives VkC„ 
in favor of VkM, similar to the case of the velocity for- 
mula in Eq. (|2.8|) . by studying the gradient of Eq. (j2.4[) 
with respect to k. The details are given in Appendix 
[Cl Here we merely record the result which facilitates the 
numerical evaluation Vt^ KR {k) in Eq. ([3TT5]) : 



f2 



ffKK 



(k) - -7m{ £ 



OJnC™(A n -A m ) 



E 

where 
V k M = 



[C t C m - (Cj C„ ) (gt C m )] (D I V k M C„ ) ' x J I V k MC„ 



(3.14) 



3G(£;k) 



At(£)| 



-v„(k) 



8G(£;k) 
d£ 



, 3Af(g) 



It is reassuring to note that for a Hermitian KKR matrix 
M, for which A = A and CjCj = Eq. ([3TT4]) reduces 
to its conventional formi 



fi™(k) = -7m{ £ e " Vk ^f^7 kjtg " } 



(3-16) 

From the point of view of the present paper, Eqs. (I3.14[) 
and (|3 . 15[) together are one of our two central formal 
results. It expresses the contribution n^ KR (k) to the 
Berry curvature in terms of the left and right eigenvec- 
tors, the group velocity, and k and £ derivatives of the 
KKR matrix. As will be demonstrated, these relations 
provide an efficient way of calculating Ct^ KR (k) similarly 
to the manner of Eq. (|2.8[) done for the group velocity. 
The main difference is that now we need the total deriva- 
tive VkM in Eq. (|3.15p instead of the partial derivative 
dM/dk used in Eq. (j2~8)) . 

Let us consider the term H^(k) = Vk X A v n (k). A de- 
tailed derivation performed in Appendix [C] finally gives 
the following expression 

n^(k) = 2v„(k)x 

r m r v \-Cl&c m -ct Ac n (clc m )]Dl-v k MC n 1 (3.17) 

1 ^ £4c m (A„-A m ) ' 



where the matrix A is defined by Eq. (|3 . 5[) . Due to the 
cross vector product, this term does not contribute to 
the Fermi surface integrals needed in calculations which 
follow Haldane's proposal^ 



Returning to the contribution fi^(k) defined in 
Eq. (|3.10p . on taking the curl, one finds 

^;(k)=/[(V k |*„k(r)| 2 )xr]dr = 

= 2JZe{/[*t k (r)(V k * nk (r)) x r]dr} . 

UJ 

(3.18) 

Then, using the KKR expansion of Eq. (|2.ip it follows 
( Appendix ICj) that 

n;(k) - 2 • v„(k) x Re{&J £ C n }- 

(3.19) 

Here the vectorial matrix r is given by Eq. (|3.8p . and the 
vectorial matrix fg is defined as 



(i £ ) QQ ,(£)=J& Q (£;r)r 



8£ 



dr . 



(3.20) 



To summarize the above discussion, we note that the 
formula for the full curvature (Eq. (|3.10[| ) together with 
Eqs. (|3.14l) . (|3.17p . and (|3 . 19[) constitutes a basis for cal- 
culating the conventionally defined fi„ (k) . Here we used 
the eigenvectors and eigenvalues of the KKR matrix and 
the matrix elements with respect to the local scattering 
states $q(£;t) given by Eqs. ([g^]) . fl3J§, and (|3~20l) . 
Clearly, these quantities are readily available in a KKR 
calculation which is aimed at computing the wavefunc- 
tions as well as the dispersion relation^ 

An important point to mention is that the projection 
of H n (k) along the group velocity 

n n (k) • v„(k) = -v n (k) • Im{ £ 



Dj,Cm(An-Aro) 



C!D k 



E 

-2v„(k) • Re{ £ 



[ClfC m -ClrC rl (ClC m )]xD 1 m (afI/dk)C n 



z5Lc m (A„ 



-A„ 



(3.21) 

does not have any terms connected with the energy 
derivative. This is then the second principle result of 
the current section. Its significance is that no Fermi sur- 
face integral contains energy derivatives. Therefore, for 
calculating the anomalous Hall conductivity according to 
Haldane's approach^, one does not need a numerical dif- 
ferentiation at all, since the partial derivative dM/dk has 
to be taken analytically according to Eq. (|2.9p . 

Up to now we discussed the conventional Abelian case 
when there is no degeneracy of the electronic states. In 
the next section we consider the Berry curvature in a 
general non- Abelian case. 



C. KKR formula for non- Abelian Berry curvature 

As it was discussed in detail by Shindou and Imura 
(Ref.El), the presence of degenerate Bloch bands makes 
the Berry curvature non- Abelian. Since two covariant 
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derivatives (with respect to k) along different axes do 
not commute with each other in the subspace spanned 
by the degenerate bands, the Abelian description fails. 
Namely, in the case of an i-fold degeneracy the Berry 
curvature is not any more a vector, but a vector-valued 
matrix in L-dimensional space labeled S. Those elements 
can be written a a 19 ' 29 

fiy(k) = i(Vk«ik| x |V k itjk}- 

-i E (V k Wik|wik) x («ik|Vk%k) , (3.22) 

I 



where indices i and j mean any two states from the set 
£ = {1, 2, L} of the degenerate states. Below we de- 
rive detailed expressions for the non- Abelian Berry cur- 
vature given by Eq. (|3.22p within the KKR method. 

For typographical simplicity, in this subsection we use 
the inner products which will always mean an integra- 
tion over the unit-cell. Then, using the Bloch theorem, 
Eq. p.22j) can be rewritten as (omitting index k for the 
wave functions) 



fiy(k) = ^(V k * ^ | x |V k tt 3 -) + (V k * 4 x r|*,) - (^|r x V k * y ) - £ {i(V k *i|*,) x (*,|V k * ,■>- 

zes (3.23) 
-(*i|r|*,) x (^ilVk^ + ^k**!^) x (^|r|^-)+i(^|r|^) x (*i|r|tf.,>} • 

Similar to the previous subsection, we can generalize our separation of the Berry curvature into the following contri- 
butions 

tin (k) = n% (k) + % (k) = n™(k) + n» (k) + % (k) . (3.24) 

Here O^- (k) splits into 

ij (k)_ \k*h CID^-K)^^-^ (3-25) 



and 



nv.(k) = -i { Vi x e Pl ie "-^'^fi ?ic ')^'^)Pl,v^Ci 



»0E SXc m (A,-A m ) 



" J ^ cJLo m (Ar-A* ) j+ ( 3 - 26 ) 

l<p. 



[Vi x Vj ] (cjA^ - E (gtAtcOCcJM-)) 



with the new matrix Af defined as 



(,A f jQQ/(tJ - J — ^ dr . (3.27) 



The last term in Eq. (|3.24j) can be written as 



n£ (k) = v, x [CjflCj - e (cjAtqxctfCi)] + 



zes _ zes 

y- [g^rgj-SiEsCglgOg/rgjlxP^VkMCi)' _ I 3 - 28 ) 

[g| -Eigs <p| rCj (C/ C m )] x V k J\7C 3 

5, t „C m (A J -A m ) ' 



A detailed derivation of these formulas is given in Ap- 
pendix |D] To get the expressions for the Abelian case, 
obtained in the previous section, one needs to consider 
the diagonal element iljj(k) and restrict the sum over I 
just to the term I = i. 



IV. RESULTS FOR THE BERRY CURVATURE 

Here we present the results for the Berry curvature at 
the Fermi surface of Al, Cu, Au, and Pt bulk crystals. All 
of them are non-magnetic materials with space inversion 
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symmetry. As was mentioned in Section I, in such a 
case the electron states are two-fold degenerate at each 
k point. In other words, they form a Kramers doublet. 
Therefore, according to the discussion of Section III, the 
non-Abelian Berry curvature is a vector- valued matrix in 
the two-dimensional space of the two degenerate bands. 

In general, each matrix element of f2(k) is gauge de- 
pendent. It would be meaningless to visualize the el- 
ements for an arbitrary gauge. A gauge-independent 
quantity is the vector TV[fi(k)], but it vanishes for 
Kramers-degenerate bands. Other gauge independent 
quantities are Tr [5 M (k)f2 M (k)] with \i = x,y,z using 
the spin matrices 5^(k) = (\&j|/3cr M |\I>j) in the subspace 
spanned by the two degenerate bands. However, these 
quantities combine already two effects stemming from the 
Berry curvature as well as the spin mixing of the wave 
functions.— Hence, some features of the Berry curvature 
may be hidden. 

Since there is no convenient gauge invariant quantity to 
plot we have chosen a physically appealing gauge. Simi- 
larly to what was discussed in Ref. 211, it is a special linear 
combination of the degenerate states such that the off- 
diagonal matrix elements of the spin operator S z = f3a z 
in the subspace of the two degenerate bands are zero. 
Such a transformation can always be performed. For 
simplicity we present the Berry curvature for one of the 
degenerate bands only, since for the diagonal elements 
the relation fin = — fi22 holds. The off-diagonal terms 
are more complicated being not purely real, but complex 
numbers. 

In Fig. [2 we compare the three separate parts con- 
tributing to the Berry curvature from Eqs. (|3.25[) . (|3.28[) . 
and ([3~26]) . The first one (Fig. [2 (a)) is the KKR 
part Q, KKR {k) that clearly has the dominant contribu- 
tion. The maximum value of the contribution fi r (k) in 
Fig. [2](b) is less than 4% of n KKR (k), and fl v (k) shown 
in Fig. [2 (c) contributes less than 2%. The same holds 
for all the other considered systems for which only the 
total Berry curvature fi(k) is summarized in Fig. m 

Here, the interesting result is that the maximum value 
for the length of the Berry curvature is largest for Al 
which is actually the lightest element with the weakest 
atomic spin-orbit coupling. However, the region of such 
a large contribution is very small and connected to points 
where the Fermi surface touches the Brillouin zone (BZ) 
boundaries. A similar effect is known for the spin-mixing 
of Bloch states on the Fermi surface of AL 21 i 30 The en- 
hancement of the spin-mixing is induced by two mecha- 
nisms. Firstly, an avoided crossing of two bands appears 
at these points. Secondly, this avoided crossing occurs 
near the BZ boundary where the spin-orbit interaction is 
already increased due to the multiband character of the 
Fermi surface of Ali 3 ^ The same explanation, connected 
with the strength of the k-dependent spin-orbit coupling, 
holds for the enhancement of the Berry curvature. Ex- 
cept for these special points, the values for the Berry 
curvature in Al are orders of magnitude smaller than for 
all the other considered metals. In addition, the avoided 




0.0 80 0.0 2.7 0.0 1.1 



FIG. 2: The length of the diagonal components of the Berry 
curvature on the Fermi surface of Au in a.u.: a) the KKR part 
fi|f (k) according to Eq. (|3.25[) . b) the contribution from 
fiii(k) given by Eq. (I3.28[) . and c) the part f2^(k) introduced 
by the energy dependence of the basis functions according to 
Eq. (EMI. 

crossings can explain the larger contributions in Pt in 
comparison to Au which is in fact heavier, but has only 
one degenerate band at the Fermi surface. Finally, Cu, 
also with one band only, has quite small contributions 
since it is a relatively light material. 

V. SPIN HALL CONDUCTIVITY 

As a first application of the Berry curvature, calcu- 
lated within the KKR formalism, the intrinsic spin Hall 
conductivity (SHC) will be presented. This quantity was 
already calculated using a Kubo formula like expression 
for the SHE&£ ' 31 i 33 and our purpose here is to validate 
our approach to the problem. 

It might be preferable for such a comparison to cal- 
culate the anomalous Hall conductivity (AHC). As it is 
known, for this quantity the Kubo-like formula and the 
semiclassical expression are formally equivalent^ How- 
ever, for nonmagnetic systems the AHC vanishes. This 
leaves us with no choice but to calculate the SHC in spite 
of two conceptual difficulties. The first of these is the lack 
of a proper definition of the spin current operator i 32 i 34 
The second one is that, even with the frequently used 
choice of the spin current operato r 6 ' 7 ' 35 , the Kubo for- 
mula for the SHC is not equivalent to the simplified semi- 
classical theory used hereM- 

In general, the AHC can be written in terms of the 
Berry curvature as 3 - ' 10 ' 37- — 

where the distribution function f n (Ep,'k) restricts the 
integral to the states below the Fermi energy Bp. 

For the SHE this formula has to be modified to account 
for the fact that a spin and not a charge current is flowing. 
In addition, the non-Abelian nature of the Berry curva- 
ture has to be taken into account.— Let us start with the 
heuristic spin-current operator j = l/2(Sv + vS)fS 
Following the simplest interpretation of the semiclassical 
wave packet dynamic o 39 ' 40 with v — > r c = — efl x E and 
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FIG. 3: The absolute value (in a.u.) of the diagonal component of the Berry curvature for the Fermi surface of several metals. 
From left to right: Al (3rd and 5th band), Cu (11th band), and Pt (7th, 9th, and 11th band). For Al we used a logarithmic 
scale to visualize the important regions. 



S — > S, we consider the anomalous velocity induced by 
an applied electric field E. Now, we have both S and fl 
as 2 x 2 matrices in the subspace spanned by the Kramers 
doublet. If we assume the electric field to be in x direc- 
tion and restrict the discussion to the spin polarization 
in z direction, then the SHC is given by 19 



to 



2/i 



E 



— 2fn(E F ,k)Tr[p n (k)S n (k)n* n (k)] . 

BZ ( 27T ) 

Here p n (k) is the density matrix which describes the wave 
packet constructed from the two degenerate states corre- 
sponding to the wave vector k and band n. As mentioned 
above, this expression is not equivalent to the Kubo for- 
mula of Refs. @, |U and [33l. The difference is induced 
by neglecting the band off-diagonal terms stemming from 
the spin operator^ Here we mean the other bands which 
are out of the considered Kramers doublet but may be 
energetically close to it. However, the Kramers dou- 
blet is treated correctly in terms of a non-Abelian Berry 
curvature^ We leave a possible influence of these simpli- 
fications to be investigated elsewhere. Here we only show 
that within such approximations one can reproduce the 
results obtained in the more rigorous approach of the 
Kubo like formula>&~— £2- To aid the emergence of physi- 
cal insight into the content of our calculations we made 
a further simplification by assuming the spin expectation 
value for the degenerate bands to be Sy(k) = ±1. This is 
equivalent to a two current model where the spin current 
is given by I s = I + — I~ . Here "+" and "— " denotes the 
current provided by \f f+ and vp - states with a positive 
or negative spin polarization, respectively 21 Thus, the 
matrix element £2* xi in Eq. (j5.2[) corresponds to 
For an incoherent superposition of two wave packets cor- 
responding to the degenerate states of the nth band the 

density matrix takes the form p n (k) = ( q ^ ~\ . There- 
fore, the SHC can be written as o z = cr+ — a~ leading 



e 



E / J^sfn(E F ,m Z n + (V (5-3) 



h ^ J (2?r) 

BZ 
E F 



h{2n) 



d£n z {£) , 



where 



n*(£) = £ / 



IS(£) 



d 2 k 

KM 



nj+(k) . (5.4) 



Here we exploited the fact that for the Kramers pair the 
condition ft£+(k) = -fi*' - (k) holds. In fact, Eq. jS3J is 
nothing else but the formula for the AHC applied for the 
"+" subband only. It is written in terms of the energy- 
resolved Berry curvature £l z (£ ) via an isosurface (IS) in- 
tegral. 

In Fig. 2] we show the SHC as a function of Ep for 
Au and Pt calculated by Eq. (I5.3[) . It is in reasonable 
agreement with the results obtained by Guo et al.— i2S 
using a Kubo formula approach. All main features in the 
energy dependence of the conductivity are reproduced. 
The conductivities at Ep are given by 470 (ficm) -1 and 
2500 (ficm) -1 for Au and Pt, respectively. 

As it is well known, the integration of the Berry curva- 
ture over the Brillouin zone is a computationally very de- 
manding task^^ This stems from the fact that fi(k) is a 
very spiky function in the crystal momentum space. Es- 
pecially, for light elements the Berry curvature turns out 
to be small everywhere except for small regions around 
avoided crossings. The reason for that is already clear 
from the article of M. Berryi He expressed the curva- 
ture of a certain band as a sum over all the other bands 
where the difference of the band energies appear in the 
denominator. The same situation occurs in Eq. (|3.25l) . 
where the eigenvalues of the KKR matrix play the role of 
the band energies. Taking this into account, it is evident 
that the Berry curvature becomes larger if two bands are 
coming close to each other. This is exactly what happens 
at avoided crossings of any kind. As was pointed out by 
Mikitik and Sharlai in Ref . lU the Berry curvature in the 
nonrelativistic case vanishes everywhere except for de- 
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FIG. 4: The energy resolved spin Hall conductivity for Au 
(top) and Pt (bottom) according to Eq. (|5,3[l . 



generacies of points or lines. In the vicinity of such degen- 
eracies the Berry curvature is a ^-distribution function. 
Adding spin-orbit coupling to the system leads, normally, 
to avoided crossings at the degeneracies, but they still 
give rise to a Berry curvature. It can be viewed as smear- 
ing out the S distribution. Importantly, the smearing is 
proportional to the strength of the spin-orbit coupling. 
It means that for light elements with a weak atomic spin- 
orbit coupling the Berry curvature is very close to the S 
function. That makes the integration quite demanding. 
This leads to the somewhat surprising situation: systems 
with stronger spin-orbit coupling and more pronounced 
effects induced by the Berry curvature can be handled 
numerically easier than systems with tiny splitting of the 
bands. 

To highlight once more the discussion above, in Fig. 
[5] the energy-resolved Berry curvature according to Eq. 
(|5.4j) is shown for Au. Clearly, even with respect to en- 
ergy Q z (£) is a very spiky function. That requires to 
use a very dense £ and k mesh as discussed by sev- 
eral authors^^ Here we used comparable numbers of 
k points to converge the Berry curvature integrals. Ac- 
tually, in Fig. [S] only two of the three contributions, 
according to the separation given by Eqs. (|3.24l) - (|3.28l) . 
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FIG. 5: The energy resolved Berry curvature for Au. The 
red (solid) and blue (dashed) curves show the separate con- 
tributions from the KKR part R ' z {£) and the dipole part 
Q\'^(£) of the Berry curvature, respectively. 



are shown. One can see that the KKR part n* KR ' z (S) 
(red solid line) dominates, whereas the part &['*(£) (blue 
dashed lines) is negligible. We should mention that the 
part Vl\i{£ ) is even much smaller and was skipped. This 
is a consequence of the above discussion related to Fig. 
[2] Thus, only the most stable part of the Berry curva- 
ture, including no numerical derivative, contributes sig- 
nificantly to the SHC. 



VI. CONCLUSION 

We have developed an efficient method to calculate the 
Berry curvature within the KKR approach applied to the 
electronic structure of solids. An unconventional scheme 
that requires to deal with a non-Hermitian KKR matrix 
is compensated by an elegant analytical differentiation of 
this matrix with respect to the crystal momentum vec- 
tor. This advantage is a feature of the screened version 
of the KKR method and should also be useful for all 
tight-binding like computational methods. The formal 
arguments starting with the local expansion of the Bloch 
function in Eq. (|2.1I) and leading to the computable for- 
mulas (|3.25l) - (|3.28[) can be readily adopted for calcula- 
tions based on other multiple scattering approaches. In 
particular, for the LMTO method the situation will be 
simplified due to the lack of any energy dependence for 
the basis functions. The efficiency and stability of the 
proposed computational procedure is shown by calculat- 
ing the Berry curvature for Al, Cu, Au, and Pt bulk 
crystals and the spin Hall conductivity for Au and Pt. 



Appendix A: Derivation of the group velocity 

Following Shilkova and Shirokovskii— we note that the 
eigenvalues of the KKR-matrix obey A n (£ n (k);k) = 
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and hence for the total derivative V k A„(£; k) we have 
' !V -- k - -V k 4(k)M|^ w= :() . (A1) 



9k 's=E n (k) 

Therefore, 



v„(k) = V k £„(k) = _^k I/ aA^) _ (A2) 



with respect to k we obtain 

dA„(£;k) _ Dl{dM/dk)C n 

Using this formula in Eq. (|A2[) one derives Eq. 



(A4) 



The next and central move is to calculate 
d\ n (£ ; k)/9k. For the case of Hermitian KKR matrix 
M(£;k) it was done in Ref. H^. Here we generalize the 
procedure to the case of a non-Hermitian matrix. With 
the definition of the left eigenvectors^ D^M = \ n Dl it Using the KKR expansion of Eq. (EH]), we obtain 
is evident that 



Appendix B: Connection .A^(k) 
To deal with Eq. (|3.2[) , we need to calculate Vk$ n k(r). 



■ ggo(r 



(A3) 



V k * tik (r) = E[V k C3(k)$ Q (r) + C3(k)v„(k)^ 
Q 



Taking the partial derivative of both sides of Eq. (|A3[) that gives us 



(Bl) 



(*nk|V k * nk ) = £ C3*(k)V k C3(k) + v n (k) £ |C3(k)| 2 / $t ; r)^%P>dr 



(B2) 



Here we have used the fact that J w &q (£; r) d $< gg £ ' r ^ dv oc 5qqi by the properties of our KKR-basis set given by 
Eq. (|2~2|) . Then for -A£(k) we can write 

^(k)=^™(k)+ l v„(k)E|C3(k)| 2 /^( f;r )££^il) d r , where A« KR (k) = i £ (k) V k Cg(k) . (B3) 
Rewriting these expressions in a matrix form, we get Eqs. (3.3)-(3.6). 



Appendix C: Abelian curvature 



Let us derive first Eq. (|3.14p for n^ KR (k) starting from Eq. (|3.13|) . Using the completeness relation given by 
Eq. (|3.12p we can perform the following expansion 



f~1 V-> Dm Ct. f-1 V-> Cn r) 

_ ^ Dm ~ C m D m 



With the Hermitian conjugate of this expansion we have 



Then, 



VkCic„xDty k c„ 



VkCtc m X-D^V k C„ _ v k c* [c m -c n (clc m )]xD f m -y k c n 

~f D 1 m C m ~~ D m C m 



(CI) 



(C2) 



(C3) 



where we have used that V k C^C n x C^V k C n vanishes since V k C„ is purely imaginary. Using Eq. (|C1[) we can 
rewrite the part of the numerator in Eq. (|C3[) as 



fc^n 



[clc m -(clc n )(cjc m )]p k 



ClD k 



In addition, due to D'-Ci = for j ^ i, we have 



D\v k Ci = -\ — 



(C4) 



(C5) 
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Therefore, finally we can write 

V VkClCrnXD^VkCn _ v 1 v [cj C m - (Cj C n ) {Cj C m )] (D I V k MCn ) * xDlV k MC n 

m mm m^n v ' k=jin k 

and end up with Eq. (|3.14p . 

For a derivation of J~J^(k) we need to take the curl of the second term in r.h.s. of Eq. (|B3|) . Taking into account 
that Vk x v„(k) = Vk x V k £„(k) = 0, one can write 

n*(k) = 2v n (k) x / m { S [C3*(k)V k ^(k)]/$ Q (£;r)^|^dr} = 2v n (k) x 7m{(7t AV k OJ . ( C 7) 

Q w 

Then, using again the completeness relation of Eq. (|3.12l) together with Eqs. (|C2[) and (|C5|) we obtain 



C\ AV k C„ = E = C-t A(7 n (7t VkC„ + E 



C} l AC m Dl 1 'V k MC n _ (ClAC n )Dl,'V k MC n (ClC m ) 
DLO m (\ n -\ m ) D 1 m C m (X n -X m ) 



Here the term C^AC n C^ V k C„ is purely real since V k C„ and C^AC„ both are purely imaginary quantities (latter 
one due to the normalization J <E>q | 2 c?r = 1). Hence we end up with Eq. (|3.17p . 

Let us consider now Eq. (|3.18[) and use the KKR expansion given by Eq. (|2.ip . Then, 

n;(k) = 2Re{ E J[C^^ Q V k (C^ Q/ ) x r]dr} = 2Re{ E J>„(k) x C^ Q r^^C^]dr}- 

QQ'uj ' QQ'ui 

(C9) 

-2i?e{ E J[C£*^r$ Q , x V k <73,]dr} = 2i?e{v n (k) x CtrgOj - 2Re{Clf x V k C„} . 

Using again the completeness relation of Eq. (|3. 121) together with Eqs. (|C2[) and (|C5p , for the second term of Eq. (|C9[) 
we can write 

r^ts v V7 r> _ c,i?c m x5^v k c„ _ v r>tY7 r> g^gr,xJ^v k jgg„(gtg m ) gtfg m x£>^v k X/c„ 
C„r x V k C„ - E qTcZ - ^„i-C„ x C n V k C„ - E Dtc m (x n -x m ) + ^ D*,c m (A„-A m ) ' 

(cio) 

Here C\yC n x C^VkCn does not contribute to Eq. (IC9|) since the quantity C|rC n is purely real while C|V k C„ is 
purely imaginary. Thus, finally we obtain Eq. (13. 19[) . 

Appendix D: Non-Abelian curvature 

We start with fi*- (k) part of the representation for the non-Abelian curvature given by Eq. (|3.24p . The first term 
contributing to this part is 

ft™(k) = iV k G\ x V k Q - i E VkCtQ x C/V k Q =zg v k ctg x^v.c^ 

(Dl) 

[V k CjC m -£ ies V k C, t C,(C, t C m )]xDt,V k MC j 



i E V k CjCi x 

ZGE 



Vk ^' = * £ m:c mi x 3 -x m) 



where we have used Eqs. (|Cip and (jC5p . According to Eq. (JUT}, we can rewrite the term in the square brackets as 

v k cj{c m - e Q(clc m )} = v k c/{E fff^fc - E (c^m)^ + e §§-5 k }} = 



(D2) 

Therefore, we end up with Eq. (|3.25p . Now we consider the second term contributing to fiy(k). Namely, 
nv.(k) = z[v, x ClAtVkQ - V j x V k CtAQ] - i E {viClAtQ x C/VtCj - Vj x VkCjQCC/AQ)} + 



* [v< x Vj ] (ctAgc^ - E (<A+q)(c+ Ac,-)} 



(D3) 
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where the matrix is defined by Eq. (|3.27|) . Here, due to Eqs. (j3. 12[) and (jCll) . we have 

v 4 x CjATVkCj = v, x E ^T7^ 1 . v j x VkCTAQ = v 3 xE J — 



'/ft /ft 

E v.-C+Attf, x C/V k C, = v 4 x E C?At<7,[ A + E if^^VkQ , (D4) 
E v, x V V G+Q(C? AC,-) = Vj x E (C/AC^VkC^ a + ^ J^Lfy . 

Hence we end up with Eq. p.26[) . Let us consider now 

n£(k) = (V k 4i x r|* 3 -) - (*i|r x V k ^> + 

+ E {(*ik|*i> x (^|V k ^-) - (V k *i|*,> x - »<*i|r|*i> x (*i|r|* 3 -)} . 

Here, due to the completeness relation given by Eq. (|3.12p . 

(V k * 4 x r|^-) - v, x CjhCj + V k C| x fCj = v t x C}f £ C 3 + E 



ts^/^. _i_ Y7. /^t v *n — -.r. v r^ts^r 1 . _i_ ^kC} D m xCl l rC j 

CLD m 



(¥,|r x V k *i> = -v, x C/ffC,- + C/r x V k C, = - Vj x + E ^"ffi v ^ 



(D5) 



(D6) 



In addition, taking into account Eq. (IC1[) . we have 

E (*i|r|*,) x C*i|Vk%) = E C|fC x [C/VkC,- + Vj-C/AC,-] = 

= EIClrQ x [JL + E ^^VkCj-vj x (CjiCiXCjACj)} , 
E (V k * 4 |* ( ) x <*,|r|* 3 -> = E [V k C|Q + Vi C|AtC] x cjfCj = (D7) 

= E{V k C|[A + E #5 m ] x C/fC,- + v, x (C/fQ)(CjAtC ; )} , 
E (*i|r|*,) x - E G\fd x C/fQ . 



Substituting these expressions in Eq. (|D5[) . we end up by the Deutsche Forschungsgemeinschaft (SFB 762). We 
with Eq. (|3.28j) . thank Sergey Ostanin who had drown our attention to 

the papers of N. A. Shilkova and V. P. Shirokovskiii 23 i 25 
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